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ABSTRACT 

The third catalog of high-energy gamma-ray sources detected by the 
EGRET telescope on the Compton Gamma Ray Observatory includes data 
from 1991 April 22 to 1995 October 3 (Cycles 1, 2, 3, and 4 of the mission). In 
addition to including more data than the second EGRET catalog (Thompson 
et al. 1995) and its supplement (Thompson et al. 1996), this catalog uses 
completely reprocessed data (to correct a number of mostly minimal errors 
and problems). The 271 sources (E >100 MeV) in the catalog include the 
single 1991 solar flare bright enough to be detected as a source, the Large 
Magellanic Cloud, five pulsars, one probable radio galaxy detection (Cen A), 
and 66 high-confidence identifications of blazars (BL Lac objects, flat-spectrum 
radio quasars, or unidentified flat-spectrum radio sources). In addition, 27 
lower-confidence potential blazar identifications are noted. Finally, the catalog 
contains 170 sources not yet identified firmly with known objects, although 
potential identifications have been suggested for a number of those. A figure 
is presented that gives approximate upper limits for gamma-ray sources at any 
point in the sky, as well as information about sources listed in the second catalog 
and its supplement which do not appear in this catalog. 

Subject headings: gamma rays: general 
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1. Introduction 

EGRET (Energetic Gamma Ray Experiment Telescope) is the high-energy gamma-ray 
telescope on the Compton Gamma Ray Observatory (CGRO). Descriptions and capabilities 
of the instrument are given by Hughes et al. (1980), Kanbach et al. (1988), Kanbach et 
al. (1989), Thompson et al. (1993a), and Esposito et al. (1998). The telescope covers 
the energy range from 30 MeV to over 20 GeV. EGRET records gamma-ray photons 
individually as electron-positron pair production events, which are processed automatically 
(with manual verification) to provide the arrival direction and energy of each photon. The 
point spread function (PSF) is energy-dependent, having a FWHM of approximately 6° 
at 100 MeV and smaller values at higher energies. The arrival time of each gamma ray is 
recorded in Universal Coordinated Time (UTC) to an accuracy of better than 100 /j,s. The 
field of view of EQRET extends to more than 30° from the instrument axis, although the 
sensitivity at angles beyond 30° is less than 15% of the on-axis sensitivity. Because of the 
low flux level of the high energy gamma rays, observing periods are typically 2-3 weeks. 

The first EGRET catalog (Fichtel et al. 1994) covered Cycle 1 of the Compton 
Gamma Ray Observatory program, which began shortly after launch of the observatory 
on 1991 April 5 and ended on 1992 November 17. That catalog included one solar flare, 
pulsars, gamma-ray bursts, a normal galaxy, active galaxies, and unidentified sources. 

In addition to source detections, it listed upper limits for other objects thought to be of 
interest. The second EGRET catalog (Thompson et al. 1995; hereafter 2EG) included, in 
addition to Cycle 1, Cycle 2, which ran from the end of Cycle 1 until 1993 September 7, and 
all of the results were reanalyzed in as consistent a manner as possible using an improved 
calculation of the diffuse radiation (Bertsch et al. 1993b; Hunter et al. 1997; Sreekumar 
et al. 1998a). A supplement (Thompson et al. 1996; hereafter 2EGS) to the 2EG catalog 
covered new sources detected during CGRO Cycle 3 or only in the sum of Cycles 1, 2, and 
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3. The approach was similar to that of 2EG, but it did not provide reanalysis of all of 
the 2EG sources. Lamb & Macomb (1997) have presented a catalog of EGRET sources 
detected above 1 GeV; some of the sources found by them are not in this catalog because 
they are below the >100 MeV significance threshold. 

This catalog is based upon a complete reprocessing of the data used for 2EG and 
2EGS, plus the incorporation of additional data from Cycle 4. As were the earlier EGRET 
catalogs, this one is for sources with E > 100 MeV, although positions have been refined 
using higher-energy subsets of the data. 

Section 2 of this paper begins with a brief summary of the observations and the 
approach to the analysis, with references to relevant papers for details. This is followed 
in section 3 with a description of the current catalog. It is similar in arrangement to 
2EG in that the sources are listed in one table (Table 4) in order of increasing right 
ascension regardless of the type of source. In Section 4 there is a table showing the status 
of sources which appeared in 2EG and 2EGS but do not appear here. In Section 5, a figure 
is presented, similar to that in 2EG, for determining an approximate upper limit for an 
arbitrary sky position. 


2. Observations and Analysis 

For scheduling purposes, the Compton Gamma Ray Observatory mission has been 
divided into cycles corresponding to proposal cycles. Cycle 1 lasted from 1991 April, 
through 1992 November, and constituted the first complete survey of the high-energy 
gamma-ray sky. Cycle 2 covered the time interval 1992 November 17 - 1993 September 7. 
Cycle 3, 1993 August 17 - 1994 October 4, overlapped slightly with the end of Cycle 2, and 
Cycle 4 extended from 1994 October 4 to 1995 October 3. Thus this catalog covers all of 
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the portion of the mission in which EGRET carried out observations with its full field of 
view. (Since the end of Cycle 4, EGRET has been operated in a narrow-field-of-view mode, 
and during a decreasing fraction of time, in order to conserve spark chamber gas lifetime.) 

Viewing periods (time intervals with spacecraft pointing fixed) have lasted from a few 
days to three weeks. Table 1 lists the start and stop dates of each viewing period and the 
number designation used to identify it. In most cases, the gaps between viewing periods are 
a few hours, the time needed to re-point and stabilize the Compton Observatory. During 
Cycle 2, several longer gaps represent the times used to prepare to re-boost the observatory 
to a higher orbit. Although there has been some coverage of all parts of the sky, the 
coverage has been far from uniform. Figure 1 displays the EGRET exposure for the sum of 
Cycles 1, 2, 3, and 4 in Galactic coordinates, in units of cm 2 s . 


2.1. Maximum Likelihood Application 


Because of the low photon detection rate and the extent of the PSF, statistical 
techniques are required to analyze EGRET data. The maximum likelihood approach 
(Mattox et al. 1996) is used to estimate point source flux densities, source locations, and 
background model parameters. The likelihood statistic of binned EGRET data is the 
product of the probability for each pixel 


£=n K 


a) 


where 


Pij = 


v 


n,j\ 


is the Poisson probability of observing n,j counts in pixel ij when the number of counts 
predicted by the model is The likelihood ratio test is used to determine the significance 
of point sources. The likelihood ratio test statistic is TS = 2(lnTi — In To), where InLi is 
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the log of the likelihood of the data if a point source is included in the model, and InLo is 
the log of the likelihood of the data without a point source. In the null hypothesis, TS is 
asymptotically distributed as xl for a source at a specific position. Mattox et al. (1996) 
show that the corresponding significance is y/TSa. Monte Carlo simulation and experience 
with flight data indicate that these techniques provide reliable results. 


2.2. Diffuse Model 

The diffuse gamma radiation consists of two components, one apparently extragalactic 
and assumed to be isotropic, the other a Galactic component. The latter is highly peaked 
along the Galactic plane but dominates over the isotropic component everywhere except at 
the highest Galactic latitudes. The diffuse model used for the generation of this catalog is 
the same one used for 2EG and 2EGS. It provides a good representation to the large-scale 
EGRET data; however, it is not necessarily ideal for every point in the sky. For this reason, 
the catalog analysis for a specific point in the sky (a 15° radius about the point being 
considered) allows the isotropic and Galactic diffuse components to be independently scaled 
from their nominal values (as was the case for 2EG and 2EGS). 


2.3. Source Detections 

The maximum likelihood program was used to examine each of the following > 100 MeV 
maps: 

each of 165 viewing periods; 

20 maps of sums of 2 - 6 viewing periods (defined in Table 6) with similar pointing 
directions and close in time; 
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the four summed maps of Cycle 1, Cycle 2, Cycle 3, and Cycle 4; 
the two summed maps of Cycles 1+2 and Cycles 3 + 4; 
the summed map of Cycles 1 + 2 + 3 + 4. 

Any excess above the diffuse model which showed a y/T S > 3 was initially retained 
as a candidate source. A cross-comparison of the candidate sources was made in order 
to eliminate cases in which the sairie source was assigned different positions in different 
observations. From this comparison, a list of candidate sources was derived. Excesses with 
y/TS > 4 (y/TS > 5 within 10° of the Galactic plane) were retained as likely catalog 
sources. Unlike the processing for 2EG and 2EGS, the below-threshold excesses with 
y/TS>3 were included in the source model while optimizing the positions and fluxes for 
the sources in the catalog, because it is unrealistic to assume there are no sources below 
the catalog threshold, and the simulations described by Mattox et al. (1996) indicate that 
excesses with y/TS > 3 are more likely to be real than statistical artifacts. One effect of 
including the below-threshold excesses in the model is to strengthen slightly the TS for 
catalog sources. Also, for sources just above the catalog threshold, the positions may be 
slightly affected (hopefully improved) by the inclusion of the below-threshold excesses in 
the processing. 

The summation of the low-significance excesses is estimated to yield an average sky 
intensity of about 10 -6 cm -2 s -1 sr -1 , an order of magnitude below the extragalactic 
diffuse intensity. Since these low-significance excesses are more numerous near the 
Galactic plane than at high latitudes, they do not represent a significant fraction of the 
Galactic/extragalactic diffuse intensity anywhere in the sky. 

As mentioned above in the Abstract, this catalog uses EGRET data which were 
reprocessed (Esposito et al. 1998) to improve the calculation of instrument sensitivity as 
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a function of energy, and to correct several minor errors in the calculation of sensitivity. 
These changes do not affect the statistical significance of source detections, although they 
can have a small effect on the source fluxes. 


2.4. Source Positions 

Experience with analysis of EGRET data has shown that the (statistically) most 
significant detection gives the best position determination for a source. A few exceptions 
exist. For example, for 3EG J 1200+2847, the most significant detection is in vp206, during 
which the source was about 31° off-axis. At such an axial angle, it is known that there are 
distortions that are corrected only approximately. Thus the position for 3EG J 1200+2847 
was derived from vp418, in which the source significance is slightly lower, but the object is 
only about 17° off-axis. 

For each candidate source, likelihood test statistic maps were constructed for the 
observation (single or summed) which produced the highest \/T S for E > 100 MeV. These 
maps used energy ranges > 100 MeV, 300 — 1000 MeV, and > 1000 MeV. Because the 
EGRET PSF is narrower at higher energies, the upper energy ranges can give better source 
position information, but only if adequate statistics are available at those energies. For the 
stronger sources, the likelihood test statistic maps for the two independent energy ranges, 
300 - 1000 MeV and > 1000 MeV, were added to improve the signal. (This is not the same 
as producing a single test statistic map for E > 300 MeV.) The three or four test statistic 
maps were compared, and the one which produced the smallest error contours was chosen 
to represent the source position, as long as VT S was greater than 4, a level chosen to reflect 
a substantial degree of confidence in the detection. 

Within each test statistic map, the source position can be determined in a variety of 
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ways. The location with the highest value of \/TS is the single most probable position 
for the source. In many cases, however, a range of locations with nearly the same level of 
confidence can be found. Often the highest test statistic is not at the center of this region. 
We have chosen, therefore, to present the “best” position as the TS-weighted centroid of 
the region enclosed by the 95% confidence contour. 

Confidence levels calculated by the likelihood mapping reflect only the statistical 
uncertainty in the position of a single source. Systematic uncertainties can also affect the 
position determination: 

1. Although the model of the diffuse radiation is assumed to be an fairly accurate 
representation of the background against which a source is seen, it is certainly not perfectly 
accurate on all size scales and for all directions in the sky. (This is the main reason the 
isotropic and Galactic diffuse components are allowed to vary independently in the source 
analysis.) 

2. The presence of nearby sources can change the apparent position of a source. 
EGRET cannot easily resolve sources within 1° of each other, and sources even 5° apart 
have substantial overlap of their PSFs. Sources in the catalog close to other sources 
must be considered to have larger positional uncertainties than the statistics alone would 
indicate. In many cases a below-threshold excess occurs near a catalog source. Including 
this excess changes the position of the catalog source, particularly if it is near the catalog 
threshold. As mentioned above, we have included such excesses with \/T S > 3 during the 
processing, although they are not included in the catalog. Catalog sources whose positions 
are influenced by such excesses are noted in the catalog as possibly source-confused (“C”). 
This is a somewhat subjective judgement; however, for each source S in the catalog, 
the catalog was searched for other catalog sources within 15°, and if there was one (or 
more) of comparable or greater intensity nearby, the source S is marked as confused. The 
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definition of “nearby” depended upon the strength of the other source(s): if all of them had 
substantially lower intensity than source S, they were ignored; if one, or a sum of two or 
more, of the other sources within 5° had intensities comparable to that of source S, S was 
marked as confused; and if one or more of the other sources was much brighter than source 
S, the confusion distance was expanded beyond 5°, to as much as 15° (when the Vela pulsar 
was the confusing source). 

Figure 2 shows sample maps of locations of four bright sources. We emphasize that 
these figures are likelihood test statistic contours, not intensity contours. In each case, the 
contours indicate the 50%, 68%, 95% and 99% statistical probability that a source lies 
within that contour (Mattox et al. 1996). 

The preprint version of this paper, as well as the public electronic version available 
on the Internet, contains such location maps for all of the sources in the catalog. Most 
of the sources lie at Galactic latitudes less than 45°, and their error contours are shown 
in Galactic coordinates. For the few sources which are substantially closer to the celestial 
equator than to the Galactic plane, the location figures are in celestial coordinate so as to 
minimize distortion. 

The scales of the location figures vary; brighter sources have smaller error contours. 
Also shown in some figures are the positions of objects which might be associated with 
the gamma-ray source. The only firm identifications, other than the solar flare, the Large 
Magellanic Cloud, and those pulsars which show gamma-ray pulsations, are radio-bright, 
flat-spectrum active galactic nuclei (blazars) which have been identified with EGRET 
sources (2EG, 2EGS, and references therein; von Montigny et al. 1995; Mattox et al. 1995; 
Mukherjee et al. 1997; Mattox et al. 1997a; Mattox et al. 1997b; Bloom et al. 1997a; Zook 
et al. 1997; Bloom et al. 1997b). In addition, we believe, based on spectral arguments, that 
the nearby radio galaxy Cen A is associated with one of the catalog sources (see note below 
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on 3EG J 1324-4314). 

The 95% contour has been chosen as representative of the EGRET statistical source 
uncertainties. The error radius quoted in the catalog is the angular radius of a circular 
cone (rather than elliptical, as in 2EG and 2EGS) which contains the same solid angle as 
the 95% contour. Mattox (in preparation) provides parameters for elliptical fits to the 95% 
contours. For lower-significance sources, the 95% contour can be extremely irregular, or 
even not closed. In those cases, the error circle radius listed is 1.62 x the radius of the 68% 
contour, and an asterisk follows the value in the table. The 1.62 multiplier is derived from 
eqn. 7 of Mattox et al. (1997a). 

For some sources, even ones that are fairly strong statistically, the center of the circle 
based on the 95% contour is far from the center of, or even outside, the 50% contour. We 
take this as indicating a possible extended or multiple source, and insert the mark “em” in 
the “Notes” column of the catalog tables. This is to a large extent a subjective evaluation. 

As noted above, systematic uncertainties make these contours somewhat optimistic 
estimates of the EGRET source location capability. We have devised no quantitative way 
to depict the effect of these systematics on the error contours. Table 2 summarizes some 
relevant information about sources along the Galactic plane whose positions are known 
well: the pulsars and the bright solar flare of 1991 June 11. The table shows the offset of 
the measured position from the true position and where the true position falls with respect 
to the confidence contours. The true positions often lie in the outer regions of the error 
contours, with the Vela pulsar position outside the 99% contour (although as the brightest 
source and most significant detection, the Vela pulsar has the smallest error contours). We 
have determined that this effect in bright sources can be reduced by using map bins smaller 
than the standard 0.5°. This increases the computation time greatly; since all of the most 
significant sources are identified with objects whose positions are well known, the smaller 
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bin size was not utilized for this work. 

These results indicate that the systematics do not pose a major problem for the source 
location capability even in these regions of significant diffuse emission or strong nearby 
sources. The error contours for many of the active galactic nuclei show that the location 
capability improves for regions away from the Galactic plane. Table 3 shows offsets from 
the true positions for a number of the AGN identified in the EGRET data. A tabulation of 
the strong AGN identifications indicates that only about 47% of them lie within their 68% 
contours. 


2.5. Catalog Construction 

Once the best position for each candidate source was determined, a new analysis of 
each viewing period or summed map was made, calculating the y/TS and the flux or upper 
limit at the fixed position. In each case, a simultaneous analysis of all the sources in the 
list was done, including the lower significance excesses (\/TS > 3) as described above, to 
obtain a flux or upper limit for each source in each observation. With three exceptions, 
the > 100 MeV flux was generated assuming a photon power law with a spectral index of 
2.0 (~ E~ 2 0 ) for estimation of the source PSF. For the three brightest sources (pulsars), 
the measured power law was used because of possible effects on nearby sources (Vela: 1.7; 
Geminga: 1.5; Crab: 2.1). 

Two independent determinations were made of the spectra of the catalog sources based 
on P1234 maps for the 10 standard EGRET energy ranges. The few significant differences 
were resolved, then one of the two sets of indices was entered into the catalog. 10 Many 

10 The two spectral analyses used the same software and general analysis approach, but 
were carried out by different people. Discrepancies were due to differences in details of the 
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sources in the catalog are too weak to yield reliable spectra or spectral indices; for those 
only a note to that effect occurs in the catalog. Spectral plots are included in the on-line 
version of the catalog for those sources having an index entry in the catalog. For those 
sources whose most significant detection was not P1234, better spectra and spectral indices 
could by obtained for those more-significant observations. 

The intense and highly structured diffuse emission along the Galactic plane makes 
sources within about |6| < 10° more subject to systematic uncertainties than those at higher 
latitudes. For this reason, we have adopted the same policy used in the earlier EGRET 
catalog of including in the catalog sources with at least one detection with vTS > 4 for 
|6| > 10°, but demanding at least one detection with \/TS > 5 for |6| < 10°. The exact 
choice of the TS level and the latitude for the transition from one level to the other is 
somewhat arbitrary: the influence of the Galactic diffuse radiation does not show an abrupt 
change. 

3. Catalog Description 

The entries in the catalog (Table 4) axe given in order of increasing Right Ascension 
of the EGRET position. For each source, multiple entries are given; the first entry is the 
one from which the source position was derived. In almost all cases, this is the detection 
with the highest statistical significance. Other entries in the table give the results for the 
summed maps for Cycle 1 (denoted Pi), Cycle 2 (P2), Cycle 3 (P3), Cycle 4 (P4), Cycles 1 
+ 2 (P12), Cycles 3 + 4 (P34), Cycles 1 + 2 + 3 + 4 (P1234), and for all viewing periods 
meeting the following criteria: (1) the source was within 30° of the EGRET pointing 

analysis; in particular, in confused regions, the likelihood analysis results seem to depend on 
the order in which sources (in a map containing many sources) are optimized. 
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direction; and (2) the exposure was large enough to derive a meaningful flux or upper 
limit. (There are two exceptions to the 30° cutoff: (a) the highest-significance detection of 
3EG J1200+2847 was in viewing period 206.0, in which the source was 31° from the axis, 
and (b) viewing periods 403.0, 403.5, 411.1, and 411.5 were carried out with EGRET in 
narrow field-of-view mode, so a cutoff of 19° from the pointing direction was used) Also 
presented are results for sums of 2 - 6 individual viewing periods close in time and with 
nearly identical pointing directions. ' For each observation (or sum of observations) with 
VrS > 2, the excess is presented as a flux with its uncertainty. (Note that for the 2nd 
EGRET catalog, \/T S > 3 was used.) For y/TS < 2, the result is presented as a 95% 
confidence upper limit, derived as described in section 3.2 of Mattox et al. (1996). Each 
analysis of a viewing period (or sum of viewing periods) was carried out independently; 
therefore, the sum of counts from individual viewing periods for a specific source does not 
necessarily match closely the counts from a summed map. 

To limit the length of the catalog, upper limits with little or no significance are deleted. 
First, if the maximum flux detected from a source is F, upper limits above F + AF are 
deleted. Second, where two or more short viewing periods have been added together, upper 
limits for the individual viewing periods have been deleted if they are greater than F + AF 
(or the upper limit) for the summed interval. 

The flux uncertainties and upper limits shown are statistical only. For all except the 
brightest sources or those in confused regions, the statistical uncertainty is larger than 
any systematic uncertainties. Systematic effects include the uncertainty in the calibration 
(Thompson et al. 1993a; Esposito et al. 1998) as a function of energy and angle within 
the instrument, and the gradual change in operating performance of EGRET as the spark 
chamber gas ages. A first order correction to this latter effect has been included. At this 
time, we recommend that a residual uncertainty of 10% should be attached to any flux 
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value, in addition to the statistical uncertainty. Sources in confused regions may also have 
significant uncertainties due to the overlapping PSFs. The fluxes for such sources should 
not be considered independent of each other, particularly in searching for time variability. 

The catalog is Table 4, in which the columns are: 

Name - based on the J2000 coordinates for the best position of the source, following 
the IAU naming convention (PASP 102, 1231). It should be emphasized that measured 
positions are not as precise as the name might suggest; 

RA and Dec - the J2000 coordinates measured by EGRET, in degrees; 

l and b - the Galactic coordinates measured by EGRET; 

$95 - the radius, in degrees, of the circle containing the same solid angle as the 95% 
confidence contour- (Asterisk * denotes that the value was obtained by multiplying the 
68% radius by 1.62. This was necessary in cases of unclosed or extremely irregular 95% 
contours.); 

F - For observations (or sums of observations) with \/T S > 2, the flux (E > 100 MeV) 
in 10 -8 photon cm -2 s -1 ; for \/T S < 2, the 95% confidence upper limit (in the same units). 
The flux F for most sources was determined assuming a photon spectral index of 2. If the 
spectral index 7 (see below) differs substantially from 2.0, some additional error should 
be assumed for the flux F. An approximate flux density at 400 MeV (~ 10 23 Hz), in units 
of picoJy, can be obtained by multiplying the value of F shown by 1.7. The fractional 
uncertainty in the flux density is AF/F. The approximations used in deriving the factor 1.7 
are discussed in 2EGS; 

AF - the la statistical uncertainty in the flux. Additional systematic error of roughly 
10% should be added in quadrature (important only for very strong detections); 
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7 - the photon spectral index in F(E) ~ E -7 . The l-a error of the index is in the 
second line. The spectral index is derived from the P1234 sum, which in many cases is 
not the most significant detection. In a few sources, the index and error entries are blank, 
which means that the P1234 sum does not give a sufficiently strong detection to produce a 
useful spectrum; 

Counts - the number of > 100 MeV photons represented by the flux or upper limit. 
The fractional uncertainty in the photon count is the same as that in F, AF/F; 

\/T S - the statistical significance of the > 100 MeV detection. It is approximately 
equal to the statistical a for a single measurement at a fixed position; 

VP - the viewing period of the specific catalog line. P1234 represents the sum of 
Cycles 1, 2, 3, and 4, etc.; 0.2+ is the sum of viewing periods 0.2 - 0.5, all of which covered 
roughly the same region of the sky during the verification phase of the mission. Entries 
ending in “+” are for sums of 2 - 6 viewing periods, defined in Table 6; 

ID - P=pulsar (indicates detection of pulsed gamma radiation); G=galaxy (LMC 
only); S=solor flare; A=active galactic nucleus; a=possible active galactic nucleus - This 
identification is questionable, either because the object has low radio flux density (< 300 
mJy) or because it lies outside the 95% uncertainty contour, sometimes even outside the 
99% contour (but within the position determination map). Note that this is similar to the 
notation used in 2EG and 2EGS (albeit looser than the definition of “marginal” used in 
2EG and 2EGS), but different from that used in the first EGRET catalog (Fichtel et al. 
1994), which distinguished identifications by their statistical significance; 

Other names - the source name from 2EG or 2EGS (or the first catalog, Fichtel et al. 
(1994), if not included in 2EG or 2EGS), or other gamma-ray references. In some cases, the 
source position has moved from that given in earlier references. For identified sources, the 
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name/names of the identified object is/are given in this column. Possible identifications are 
followed by 

Notes - @=see note in text; E=extended source (applies only to the Large Magellanic 
Cloud); em=possibly extended source or multiple sources (based on source location 
maps inconsistent with a single point source or poor fit to the calibrated PSF, from the 
observation or sum of observations presented in the first entry for the source); C=source 
confusion may affect flux, significance, or position. (Below-threshold excesses are considered 
in assigning this symbol, so some weak sources may be designated as confused despite 
having no other catalog sources nearby, e.g., in Figure 4.) Sources with no entry (other 
than @) in this column are consistent with the EGRET PSF for a single source; 

References - previous references to the source identification; in many cases, only the 
most recent reference(s) are shown; 

z - redshift (AGNs only). 

In many instances, redundant entries for a source are omitted from Table 4. For 
example, if only one viewing period contained a specific source during Cycle 1, PI is 
omitted from the listing for that source; all of the available information is contained in the 
listing for the viewing period. Similarly, if a specific source was in the field of view several 
times during Cycle 3, but was not observed during Cycle 4, there is no listing for the sum 
of Cycles 3 and 4 (denoted as P34 in the table). 

It is apparent that a larger fraction of the sources in this catalog are noted as being 
potentially confused than in 2EG and 2EGS. This is not only because the catalog contains 
substantially more sources than 2EG and 2EGS, but also because for weak sources we have 
considered the below-threshold excesses as confusing. 
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4. Second Catalog Sources Not in the Third Catalog 

As a result of the re-analysis of the data, a number of the sources which appeared in 
2EG and 2EGS do not appear in this catalog. These sources are scattered throughout the 
sky. In most cases, the sources from 2EG and 2EGS which were affected enough by the 
re-analysis to be dropped from the third catalog had statistical significances which changed 
from just above the catalog threshold to just below it. These “lost” sources are listed in 
Table 6, showing the \/T S in 2EG/2EGS and the maximum \/TS found for the source 
in the current analysis. One marginal AGN identification listed in 2EG (1317+520) was 
dropped in this manner. In all cases, excesses are still seen, but with significance below the 
threshold for the current catalog. 


5. Upper Limits for Sources not Detected 

Instead of providing upper limits for specific sources, upper limits have been calculated 
for a 1° x 1° grid on the sky. The result is shown in Figure 3, where the upper limit is a 
95% confidence limit in units of photons (> 100 MeV) cm -2 s _1 . For comparison with this 
figure, the faintest source in the catalog with a significance \/TS > 4 has a flux of (6. 2±1.7) 
10 -8 cm -2 s -1 . 

The sources in this catalog must be considered in the upper limit estimate. As noted 
above, EGRET cannot easily resolve sources within 1° of each other, and the PSF of 
EGRET is large enough to influence sources several degrees away. For this reason, upper 
limits near identified catalog sources may be underestimated. (On the other hand, an 
unidentified source nearby could actually be the object in question, in which case the upper 
limit from Figure 4 would be meaningless.) Therefore the regions around catalog sources 
are blanked out in Figure 4. For any object within a blanked region, it is recommended 
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that the source flux itself be taken as a conservative upper limit. 


6. Transient Sources 

3EG J 1837-0423 is a transient Galactic plane source (Tavani et al. 1997). It appears 
in this catalog because it was above the catalog threshold in Viewing Period 423.0. 

GRO J 1125-6005 is a transient source identified with Cen X-3 by Vestrand et al. 
(1997). Its VTS = 4.7 in this analysis is below the catalog threshold for a low-latitude 
source. However, discovery of pulsation at the Cen X-3 spin period with 3cr significance 
provides considerably greater confidence for this identification. 

Among the EGRET-detected blazars, there is ample evidence for transient behavior 
and extreme variability, by as much as a factor of 100 in flux (Mattox et al. 1997b). Less 
dramatic variation is seen in many of the EGRET-detected blazars, including some of the 
relatively weak ones which have been detected only once. 

Thompson et al. (1997) have searched for transient sources in the Galactic anticenter, 
but have found no strong evidence for such. Bloom et al. (1997c) have searched for 
previously unknown transient sources at high latitudes, but have also found none (although 
some previously unknown instances of variability in known sources were found). The 
search for previously unknown Galactic plane transients is continuing (Thompson, private 
communication). 


7. Notes on Individual Sources 

3EG J0010+7309 - suggested (Lamb & Macomb 1997) as possibly associated with the 
SNR CTA 1. Brazier et al. (1998) have suggested that this source might be a Geminga-like 
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pulsar associated with CTA 1. 

3EG J0222+4253 - Verbunt et al. (1997), Kuiper et al. (1998), and Hermsen et al. 
(1998) have shown indications from EGRET data for pulsed gamma-ray emission between 
100 and 1000 MeV from PSR 0218+42, which is 1° from 3C 66A, and 1° from the catalog 
position. Those authors conclude that the data are consistent with the pulsar being the 
counterpart below 1 GeV (based on timing and spatial analysis) and 3C 66A above 1 GeV 
(spatial analysis). The analysis done for this catalog (> 100 MeV) found only one source 
near this position. The catalog position given is based on the > 1 GeV map, which gives 
the smallest error contours; the contours of that map agree well with the position of the 
BL Lac object 3C 66 A. The 100 — 300 MeV map, however, seems to exclude 3C 66 A, 
and gives TS=0 at its position. The 100 — 300 MeV map is consistent with all of the 
flux in that energy range coming from the pulsar. The spectral index (2.01) given in the 
catalog assumes only one source; if there are actually two, the 3C 66A spectrum will be 
substantially harder and the PSR 0218+42 spectrum will be correspondingly softer. 

3EG J0530-3626 - considered a strong identification with 0521-365 in 2EG. However, 
the position in this catalog, based on a much stronger detection in Cycle 4, is well outside 
the 99% confidence contour. 

3EG J0542+2610 - suggested (Sturner & Dermer 1995) as possibly associated with the 
SNR G180.0-1.7 (S147). 

3EG J0617+2238 - suggested (Sturner & Dermer 1995; Esposito et al. 1996) as possibly 
associated with the SNR G189. 1+3.0 (IC 443). 

3EG J0628+1847 - Since this source is only 1.6° from the very bright Geminga pulsar, 
it could conceivably be an artifact resulting from imperfections in the PSF. 


3EG J0631+0642 - suggested (Sturner & Dermer 1995; Esposito et al. 1996) as possibly 
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associated with SNR G205.5+0.5 (Monoceros). 

3EG J0824-4610 - suggested (Sturner & Dermer 1995) as possibly associated with the 
SNR G263 9-3.3 (Vela). 

3EG J0827-4247 - close to SNR G260.4-3.4 (Puppis A), suggested (Esposito et al. 
1996) as a possible source of high energy gamma radiation. 

3EG J0834-4511 - The identification of this source, the brightest in the sky on average, 
is not in doubt, because essentially all of the gamma radiation is pulsed at the frequency of 
PSR 0833-45. The large discrepancy between the position contours and the true position of 
the pulsar is due to the 0.5° binning of the standard maps, compared with the 0.25° size 
of the position map. A similar but less dramatic effect is seen in the Geminga pulsar. This 
effect is not apparent for the Crab pulsar because its steeper spectrum makes the position 
determination less accurate. Also, the apparent large time variation of the period-averaged 
flux suggested by the individual observations is not supported by careful analysis of the 
data (Ramanamurthy et al. 1995; Kniffen, private communication). 

3EG J0841-4356 - Since this source is only 1.7° from the very bright Vela pulsar, it 
could conceivably be an artifact resulting from imperfections in the PSF. 

3EG J 1102-6103 - suggested (Sturner &; Dermer 1995) as possibly associated with the 
SNR G291. 0-0.1 (MSH 11-62). 

3EG J 1222+2841 - BL Lac object 1219+295 = W Comae = ON +231 is well outside 
the 99% confidence contour. However, the identification with this object is considered 
strong based on the position found for E > 1 GeV by Lamb Macomb (1997). 

3EG J 1324-4314 - good position agreement with Cen A, the nearest and brightest radio 
galaxy. Although this is the only candidate detection of a close radio galaxy by EGRET, 
we believe the identification is robust, since the spectrum of this source agrees well with the 
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extension of the OSSE/COMPTEL spectrum (Sreekumar et al. 1998b; Sreekumar et al. 
1999). That spectrum, in turn, connects well with the spectrum at lower X-ray energies, for 
which the Cen A identification is clear. 

3EG J 1410-6147 - suggested (Sturner & Dermer 1995) as possibly associated with the 
SNR G312.4-0.4 . 

3EG J 1627-2419 - Although the gas in Rho Ophiuchus is included in the background 
model, the scale of the variations is much smaller than the pixel size in the EGRET maps. 
This may lead to apparent sources. 

3EG J 1638-5155 - Roughly identified with 2EG J1648-5042; however, because of the 
appearance in this catalog of three nearby sources, as well as two below-threshold excesses, 
the position shown here is almost 2° away from the 2EG position. The suggestion by 
Mattox et al. (1997) of identification with radio source MRC 1646-506 = PMN J1650-5044 
appears much less convincing with the new position. 

3EG J 1800-2338 - suggested (Sturner & Dermer 1995; Esposito et al. 1996) as possibly 
associated with the SNR G006. 6-0.1 (W28). 

3EG J 1823-1314 - suggested (Sturner & Dermer 1995) as possibly associated with the 
SNR G018.8+0.3 (Kes 67). 

3EG J 1856+01 14 - suggested (Sturner & Dermer 1995; Esposito et al. 1996) as 
possibly associated with the SNR G034. 7-0.4 (W44). De Jager &; Mastichiadis (1997) have 
developed a model for the presumed high-energy gamma-ray emission from W44. 

3EG J 1903+0550 - suggested (Sturner & Dermer 1995) as possibly associated with the 
SNR G040.5-0.5 . 

3EG J2020+4017 - suggested (Sturner & Dermer 1995; Esposito et al. 1996) as possibly 



- 24 - 


associated with the SNR G078. 2+2.1 (7 Cygni). 

3EG J2035+4441 - close to SNR W63, suggested (Esposito et al. 1996) as a possible 
source of high energy gamma radiation. 


8. Discussion 

In regions of the sky where two or more sources have overlapping PSFs, it is often 
difficult to resolve the number of sources and their locations. This analysis is made 
more difficult by time-varying sources and nearby sources with different energy spectra. 
The catalog entries in such regions (i.e., the sources marked “C”) cannot, therefore, be 
considered unique solutions. (See 2EG for an example.) 

The \/TS values in the catalog represent the statistical significance for a single source 
at the given position. Because a large number of observations are included in the catalog 
(165 individual viewing periods and 27 summed maps), the number of trial positions is 
large. Following the method of Mattox et al. (1996), we estimate that the number of 
high-latitude sources (|6| > 10°) in this catalog that are statistical artifacts is between 
11 and 21. Thus high-latitude sources near the catalog threshold, especially those with 
indications of confused or extended emission, should be treated with considerable caution. 

There should be essentially no statistical artifacts among the low-latitude sources. 
Systematics may lead to some artifacts, but we do not know how to estimate the number 
of such. Presumably, sources near the \/TS = 5 threshold are suspect, but defects in the 
diffuse background map could produce apparent sources considerably more significant. 

It is known that the EGRET PSF used in this work is not perfect, first because of 
statistical limitations of the calibration data, and second, because it assumes a single power 
law for all energies above 100 MeV. (There is work in progress to derive a better PSF from 
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in-flight data, but it was not yet available for use here.) It has been shown by modeling 
(Willis, private communication) that this might lead to artifacts near extremely bright 
sources such as the three bright pulsars Crab, Geminga, and (especially) Vela. Two sources 
very near the Vela and Geminga pulsars are included in the catalog. Despite the fact that 
we have done several crosschecks that appear to confirm the reality of the nearby weaker 
sources, they should be treated with extreme skepticism. In the catalog, those sources are 
marked with a special note. 

Many of the high latitude sources, especially the blazars, are strongly time-variable. 
The pulsars show no strong time variability (Ramanamurthy et al. 1995). Some unidentified 
Galactic sources appear to be time variable (McLaughlin et al. 1996); a few (but not 
all) should be blazars seen through the Galactic disk. Other Galactic sources may be 
pulsars, either unpulsed emission from radio pulsars or radio-quiet pulsars like Geminga; 
see Thompson et al. (1994), Romani & Yadigaroglu (1995), and Kaaret & Cottam (1996) 
for discussions) or supernova remnants (Sturner & Dermer 1995; Esposito et al. 1996). 
Some of the steady high-latitude sources could be nearby pulsars (e.g. Mukherjee et al. 
1995a). As was true for the first EGRET catalog (Fichtel et al. 1994) as well as 2EG and 
2EGS, the unidentified sources show no counterparts at other wavelengths which lead to 
clear identifications. The identification of new classes of astrophysical objects capable of 
producing such high-energy nonthermal radiation remains a topic of great interest. 


9. Summary 

This third EGRET catalog contains 271 sources. Their locations are shown in Figure 5 
in Galactic coordinates; the size of the symbol represents the highest flux seen for the 
source. The 80 sources with |6| < 10° consist of five pulsars, one solar flare (bright enough 
to be detected in the maps as a source), and 74 unidentified sources. The 181 sources 
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with \b\ > 10° include the LMC, 66 high-confidence identifications with blazars, one likely 
detection of a radio galaxy (Cen A), 27 sources which may be AGN, and 96 unidentified 
sources. 

Fourteen sources reported in 2EG and 2EGS do not appear in this catalog because 
their statistical significance is below the catalog threshold in this analysis. One of these was 
marginally identified with an AGN in 2EG. 

The presence of a high fraction of unidentified sources in the catalog indicates the 
need for continuing studies, not only with EGRET and future gamma-ray missions, but 
across the electromagnetic spectrum, in order to find counterparts capable of producing this 
high-energy radiation. 


10. On-Line Catalog Information 

The catalog tables are available by anonymous ftp from 
ftp://gamma.gsfc.nasa.gov/pub/cat3/. The format is plain ASCII. The figures, 
including the full set of source location maps, are also available, in PostScript format, at 
the same location. The map of upper limits for any point in the sky is available there 
in FITS and PostScript formats. For those sources whose spectral index appears in the 
catalog, spectral plots are included in PostScript form, as well as FITS files containing the 
information in the spectral plots. 

Most or all of the information listed above is also available 
through the Web pages of the CGRO Science Support Center (COSSC), 
http://cossc.gsfc.nasa.gov/cossc/egret/egretform.html, where there are also some 
simple sorting and analysis tools. 

The list of below-threshold excesses will not be distributed publicly, but can be 
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obtained for appropriate applications by contacting the first author. 
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ll. Figure Captions 

Fig. 1 — EGRET sky exposure in units of 10 8 cm 2 s (> 100 MeV) for the sum of Compton 
Gamma Ray Observatory Cycles 1, 2, 3, and 4 (1991 April 20 - 1995 October 3). The 
contours are spaced at intervals of 2 x 10 8 . 

Fig. 2.- Location maps for four bright sources in this catalog. The contours represent 
50%, 68%, 95%, and 99% statistical probability that a single source lies within the given 
contour. A cross marks the location of objects either identified with or possibly related to 

the gamma-ray source. 

pig 3 — Map showing approximate upper limits for sources at (almost) any position in 
the sky, in units of 10" 8 photons (> 100 MeV) cm -2 s -1 , based on the summed exposure 
of Cycles 1, 2, 3, and 4, and accounting for the effect of diffuse emission as well as possible 
sources below the catalog threshold. Within 1° of catalog sources the map is black, indicating 
that this technique is not useful in those small regions. 

Fig. 4. — The third EGRET source catalog, shown in Galactic coordinates. The size of the 
symbol represents the highest intensity seen for this source by EGRET. 
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1992 

Dec 22 

77.85 

0.69 

1992 

Dec 29 

294.70 

61.88 

1993 

Jan 05 

294.46 

61.58 

1993 

Jain 12 

294.70 

61.88 

1993 

Feb 02 

314.06 

31.51 

1993 

Feb 09 

307.39 

20.75 

1993 

Feb 22 

0.24 

-34.01 

1993 

Feb 25 

355.62 

6.28 

1993 

Mar 09 

125.86 

-4.70 

1993 

Mar 23 

83.74 

11.67 

1993 

Mar 29 

182.63 

-8.22 

1993 

Apr 01 

355.62 

6.28 

1993 

Apr 06 

311.66 

22.89 

1993 

Apr 12 

140.75 

38.11 

1993 

Apr 20 

311.66 

22.89 

1993 

May 03 

151.41 

71.26 


1993 

1993 

1993 

1993 

1993 


May 06 
May 13 
May 24 
May 31 
Jun 03 


350.10 

298.09 

187.52 

157.79 

359.14 


-44.63 

-5.88 

70.63 

-0.09 


1993 

Jun 14 

298.09 - 

-44.63 

1993 

Jun 29 

355.00 

5.00 

1993 

Jul 13 

148.11 

41.22 

1993 

Jul 27 

149.86 

42.69 

1993 

Jul 30 

276.66 

-2.27 

1993 

Aug 03 

278.79 

1.44 

1993 

Aug 10 

22.22 

-13.08 

1993 

Aug 11 

5.00 

5.00 

1993 

Aug 17 

5.00 

5.00 

1993 

Aug 24 

263.59 

-2.74 

1993 

Sep 07 

347.50 

0.00 

1993 

Sep 09 

89.13 

7.82 

1993 

Sep 21 

1.41 

9.26 

1993 

Sep 22 

277.21 

12.83 

1993 

Oct 01 

89.13 

7.82 

1993 

Oct 04 

64.33 

25.27 

1993 

Oct 19 

89.13 

7.82 

1993 

Oct 25 

278.20 

66.70 

1993 

Nov 02 

277.71 

62.70 

1993 

Nov 09 

277.60 

58.70 

1993 

Nov 16 

268.69 

69.24 

1993 

Nov 19 

283.22 

74.65 

1993 

Dec 01 

283.22 

74.65 

1993 

Dec 13 

195.14 

4.27 

1993 

Dec 15 

283.70 

74.50 

1993 

Dec 20 

283.70 

74.50 

1993 

Dec 27 

280.50 

70.70 

1994 

Jan 03 

289.30 

78.70 

1994 

Jan 16 

304.18 

-0.99 

1994 

Jan 23 

304.18 

-0.99 

1994 

Feb 01 

309.52 

19.42 

1994 

Feb 08 

68.44 

-0.38 

1994 

Feb 15 

181.44 

-2.64 



TABLE 1 cont. 

EGRET Viewing Periods 


TABLE 1 cont. 

EGRET Viewing Periods 


VP 


start 


end 


1 


321.5 

317.0 

319.0 

320.0 

319.5 

323.0 

322.0 

324.0 

325.0 

326.0 

327.0 

328.0 

329.0 

331.0 

330.0 

331.5 

332.0 

333.0 

335.0 

334.0 

335.5 

336.0 

336.5 

337.0 

338.5 

339.0 

401.0 

402.0 

402.5 


1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 


Feb 15 
Feb 17 
Mar 01 
Mar 08 
Mar 15 
Mar 22 
Apr 05 
Apr 19 
Apr 26 
May 10 
May 17 
May 24 
May 31 
Jun 07 
Jun 10 
Jun 14 
Jun 18 
Jul 05 
Jul 12 
Jul 18 
Jul 25 
Aug 01 
Aug 04 
Aug 09 
Aug 31 
Sep 20 
Oct 04 
Oct 18 
Oct 25 


1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 

1994 


Feb 17 
Mar 01 
Mar 08 
Mar 15 
Mar 22 
Apr 05 
Apr 19 
Apr 26 
May 10 
May 17 
May 24 
May 31 
Jun 07 
Jun 10 
Jun 14 
Jun 18 
Jul 05 
Jul 12 
Jul 18 
Jul 25 
Aug 01 
Aug 04 
Aug 09 
Aug 29 
Sep 20 
Oct 04 
Oct 18 
Oct 25 
Nov 01 


181.44 

158.48 

143.99 

83.09 

146.42 

356.84 

197.01 

15.03 

147.04 

195.92 

82.86 

64.87 

253.40 

64.87 

18.00 

64.87 

18.00 

64.87 

253.40 

9.01 

253.40 

88.37 

340.43 

205.00 

263.59 

4.05 

113.90 

310.28 

306.73 


-2.64 
45.38 
28.02 
45.47 
26.02 
-11.29 
58.62 
5.63 
-9.04 
58.31 
-49.56 
-0.03 
142.00 
-0.03 
0.00 
-0.03 
0.00 
-0.03 
-42.00 
-8.38 
-42.00 
-46.83 
2.86 
-13.00 
-2.74 
40.40 
6.22 
-4.99 
-3.78 


VP 


start 


end 


1 


403.0 
403.5 

404.0 

405.0 

406.0 

407.0 

408.0 

409.0 

410.0 
411.1 

411.5 

412.0 

413.0 

419.1 

415.0 

418.0 

419.5 

420.0 

421.0 

422.0 

423.0 

423.5 

424.0 

425.0 

426.0 

427.0 

428.0 

429.0 


1994 

1994 

1994 

1994 

1994 

1994 

1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 


Nov 01 
Nov 09 
Nov 15 
Nov 29 
Dec 13 
Dec 20 
Jan 03 
Jan 10 
Jan 24 
Feb 14 
Feb 21 
Feb 28 
Mar 07 
Apr 04 
Apr 11 
Apr 25 
May 09 
May 23 
Jun 06 
Jun 13 
Jun 20 
Jun 30 
Jul 10 
Jul 25 
Aug 08 
Aug 22 
Sep 07 
Sep 20 


1994 

1994 

1994 

1994 

1994 

1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 
1995 


Nov 09 
Nov 15 
Nov 29 
Dec 07 
Dec 20 
Jan 03 
Jan 10 
Jan 24 
Feb 14 
Feb 21 
Feb 28 
Mar 07 
Mar 21 
Apr 11 
Apr 25 
May 09 
May 23 
Jun 06 
Jun 13 
Jun 20 
Jun 30 
Jul 10 
Jul 25 
Aug 08 
Aug 22 
Sep 07 
Sep 20 
Sep 27 


58.15 

206.81 

7.23 

306.67 

336.33 

334.33 
305.11 

274.68 
82.19 

145.10 

143.30 

185.34 
191.77 

207.35 
275.72 
158.08 
211.86 
198.21 
355.33 
355.41 

2.59 

345.74 
312.67 
137.35 
184.50 

153.75 
270.59 

18.34 


37.52 
35.82 
-73.43 
56.54 
67.23 
62.98 
57.06 
-39.17 
-32.57 
23.85 
22.72 
0.67 
-3.35 
-19.05 
-24.01 
65.84 
-17.56 
-18.26 
0.38 
-0.37 
- 0.21 
13.45 
19.04 
-47.31 
-5.87 
-9.99 
82.47 
3.98 



TABLE 2 



Solar Flare 
Crab Pulsar 
Geminga Pulsar 
Vela Pulsar 
PSR B1055-52 
PSR B1706-44 


0.39 

0.06 

0.08 

0.08 

0.17 

0.20 


68 % 

95% 

>99% 

>99% 

95% 

>99% 


TABLE 3 

Source Location Accu racy for Selected AGN 


0208-512 

0528+134 

1406-076 

1633+382 

Mkn 421 

3C273 

3C279 

3C454.3 


0.11 

0.15 

0.23 

0.12 

0.06 

0.12 

0.07 

0.13 


95% 

68 % 

95% 

68 % 

50% 

50% 

95% 

50% 
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Third EGRET Source Catalog 
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TABLE 5 

Sources from the Second EGRET Catalog and 
Supplement not appearing in the Third Catalog 


Source 

VTS in 2 EG 

present analysis 

2EG J0403+3357 

4.5 

3.2 

2EG J0426+6618 

4.5 

3.2 

2EGS J0500+5902 

4.0 

3.3 

2EGS J0552-1026 

4.3 

3.5 

2EG J 1 136 — 0414 

4.1 

3.2 

2EGS J 1236 — 0416 

4.2 

3.9+ 

• 2EG J1239+0441 

6.3 

3.9 1 

2EG J1314+5151 

4.0 

3.6 

2EG J 1430+5356 

4.1 

3.8 

2EG J 1443-6040 

5.2 

4.3 

2EG J 1631 — 2845 

6.0 

3.9 2 

2EG J1709-0350 

4.3 

3.9+ 

2EG J 1815+2950 

4.0 

3.8 

2EG J2027+1054 

4.4 

2.9 


founts split between 3EG J1236-f0457 and a source below the 
catalog threshold 

2 counts split between 3EG J 1625-2955, 3EG J 1638-2749, and a 
source below the catalog threshold 
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